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ABSTRACT. In recent years, a remarkably simple physical picture of the 
Lyman-alpha forest has emerged, which allows detailed predictions to be made 
and turns the forest into a powerful probe of cosmology. We point out ways in 
which such a picture can be tested observationally, and explore three areas in 
which the Lyman-alpha forest can yield valuable constraints: the reionization 
history, the primordial mass power spectrum and the cosmological constant 
or its variants. The possibility of combining with other high redshift observa- 
tions, such as the Lyman-break galaxy surveys, to provide consistency checks 
and complementary information is also discussed. 



1 Introduction 

The Lya forest was predicted 
and observed in the 60s (Gunn & 
Peterson 1965, Bahcall & Salpeter 
1965, Lynds & Stockton 1966, Bur- 
bidge et al. 1966, Kinman 1966; 
see Rauch 1998 for further ref.). 
Since then, there have been many 
attempts to place the study of the 
forest within the framework of cos- 
mological structure formation theo- 
ries (e.g. Doroshkevich & Shandarin 
1980, Rees 1986, Bond et al. 1988, 
McGill 1990, Bi et al. 1992). Recent 
numerical simulations lent support 
to some of these ideas, while clar- 
ifying the nature of the forest and 
allowing detailed predictions to be 
made (e.g. Cen et al. 1994, Hernquist 
et al. 1995, Zhang et al. 1995, Pe- 
titjean et al. 1995, Miralda-Escudc 
ct al. 1996, Bond & Wadsley 1997, 
Theuns et al. 1998). We explain in 
this contribution the physical picture 
of the forest that has emerged from 
these and other semi-analytical work 
(§^), discuss observational tests of 
this picture (§^|), and point out three 
examples of how such a picture al- 
lows us to "mine" valuable cosmolog- 
ical information from the forest (§«)■ 

2 Simplicity of the Forest 

Numerical 
simulations indicate that the low 
column-density (JVhi ~ 10 14 ' 5 cm — 2 ) 
Lyc? forest (for redshift zZ2), which 
occupies a large fraction of any 
given quasar absorption spectrum, 



arises from density (and velocity) 
fluctuations of a smooth intergalac- 
tic medium. Regions of enhanced 
density (in redshift-space) in this 
smoothly fluctuating medium natu- 
rally suffers stronger absorption, cre- 
ating features like absorption lines, 
which cannot necessarily be inter- 
preted as due to discrete, well- 
isolated clouds, as in traditional the- 
ories. 

The dynamics of such a medium 
is simple: gravitational instabil- 
ity on large scales, and smooth- 
ing due to finite gas-pressure on 
small scales (Reisenegger & Miralda- 
Escude 1995, Bi & Davidsen 1997, 
Hui, Gnedin & Zhang 1997). Hence, 
the dark matter and baryon com- 
ponents trace each other on large 
scales. The fluctuations giving rise to 
the low JVhi forest are mildly non- 
linear, with the overdensity Sp/pi: 5, 
where p is the baryon density and 
p its mean. The thermodynamics 
is well-understood. Shock-heating is 
minimal and a tight temperature- 
density relation exists as a result of 
three processes, recombination cool- 
ing, photoionization heating and adi- 
abatic heating/cooling: T = Top 7-1 
where To and 7 depend on reioniza- 
tion history. 

The optical depth r, which is re- 
lated to the probability of transmis- 
sion (/) by / = exp[-r], is (e.g. Hui, 
Gnedin & Zhang 1997) 

t(u ) = J Ap a (x)W[u(x)-u ]dx{l) 

where uq is the velocity of obser- 
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vation, x is the comoving distance 
along the line of sight, u is the to- 
tal (Hubble plus peculiar) velocity 
as a function of position, and W is 
the Voigt-profile given by W[A«] = 
(■Kb^)-°- 5 exp[-Au 2 /b^] with b T = 
y/ 2kgT/m p where kg is the Boltz- 
man constant and m p is the proton 
mass. The neutral hydrogen (which 
is what causes the Lya absorption) 
density is related to the baryon den- 
sity p by Ap a , where 1.6~a~1.8 
for reasonable reionization histories 
(Croft et al. 1997, Hui & Gnedin 
1997). The constant A depends on a 
number of parameters, among them 
the cosmological mean baryon den- 
sity, and the ionization background 
intensity J whose size is uncertain by 
an order of magnitude observation- 
ally. 

To predict the observational prop- 
erties of the forest, all one needs 
to know are then I. the density 
and velocity distributions of the 
gas, which can be obtained by 
hvdrodynamic simulations (ref. in 
§pJ), analytical/semi-analytical ap- 
proximations made possible by the 
mild nonlinearity of the fluctuations 
(e.g. Reisenegger & Miralda-Escude 
1995, Bi & Davidsen 1997, Gnedin 
& Hui 1996, Hui et al. 1997) or alter- 
native efficient numerical techniques 
(e.g. Gnedin & Hui 1998 |, Croft 
et al. 1998); II. the temperature- 
density relation which can be ob- 
tained using semi-analytical meth- 
ods (Hui & Gnedin 1997); III. the 
constant A, which can be fixed by 
matching the theoretical mean trans- 
mission with the observed one (Croft 
et al. 1998). 

Lastly, a few words on the higher 
column-density absorption systems: 
they are likely to arise from more 
nonlinear objects such as collapsed 
halos or galaxies. Such objects seem 
to dominate the absorption spectra 
at low redshifts (z 1) (e.g. Lanzctta 
et al. 1995), but not at higher red- 
shifts (see Rauch 1998 for a review). 



1 A regular Particle-Mesh N-body 
code can be modified to compute the 
effective potential due to pressure in 
addition to the usual gravitational 
potential, taking advantage of the 
known temperature-density relation. 



3 The Cosmic Gold Mine 

3.1 Reionization History 

Here we explore some effects of 
the reionization history on a particu- 
lar statistic called the b-distribution, 
which has been measured by several 
groups (e.g. Hu et al. 1995, Lu et al. 
1996, Kirkman & Tytler 1997). The 
b-parameter refers to the width of a 
Voigt-profile-fit of an absorption line. 
The number of lines as a function of 
width is then the b-distribution. The 
picture outlined in §U provides a sim- 
ple interpretation of the measured 
widths. According to eq. (|l]), a peak 
of p a in velocity-space which is nar- 
rower than the thermal broadening 
width fey will appear as a peak of r 
with the shape of a Voigt-profile and 
a width of fey. Hence, the measured 
b-parameter provides a direct indi- 
cation of the temperature of the gas. 
On the other hand, structure forma- 
tion theories predict and allow large- 
scale fluctuations where a peak of 
p a is much wider than the local fey, 
which means according to eq. (|j|) the 
corresponding peak in r is no longer 
given by a Voigt-profile shape and its 
width reflects more the scale of the 
fluctuation rather than the temper- 
ature of the gas i.e. the measured b 
can no longer be equated with the 
thermal broadening width. 

In more quantitative terms, we 
can always perform the following ex- 
pansion around a given absorption 
peak: 



t(u) = exp [lnr(u)] 



(2) 



r(iip) exp 



-[lnr]"(n-« p ) 2 



where u p is the velocity coordinate 
of the line center, and the prime 
denotes differentiation with respect 
to u and the second derivative is 
evaluated at u p . The first deriva- 
tive vanishes because r is at a lo- 
cal extrcmum. This expansion gives 
none other than the Voigt profile it- 
self: oc exp[— (u — u p ) 2 /b 2 ]. Hence, 
the fitted b-parameter would be fe = 
y/~ 2/[lnr]". Assuming most fitted- 
absorption lines do arise from peaks 
in r, and assuming lnr is Gaussian 
random (as in the case of linear fluc- 
tuations, or in the lognormal model), 
it can be shown that the normalized 
distribution of lines is given by (Hui 
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& Rutledge 1997) 

dN/db = 4(bi/b 5 )exp[~bi/b i ] (3) 
where 

fe4 = 2/<([lnr]") 2 > and (} refers to 
ensemble-averaging. This model nat- 
urally explains the salient features of 
the observed b-distribution: a sharp 
lower cut-off and a long high-b tail, 
which could be hard to explain in 
traditional theories of the forest (sec 
Fig. 18 of Bryan et al. 1998 for a com- 
parison of the above dN/db with ob- 
servations). In other words, high-b- 
peaks correspond to gentle fluctua- 
tions which are unsuppressed (hence 
the long tail), while low-b-peaks are 
sharp fluctuations which are statisti- 
cally rare. 

This model also tells us what de- 
termines the b-distribution, through 
the parameter b a : I. the (dimen- 
sionlcss) average amplitude of the 
fluctuations (how nonlinear the field 
is; see below), and II. the three 
smoothing scales in the problem: the 
observation resolution, the average 
thermal broadening scale and the 
baryon-smoothing-scale due to finite 
gas pressure; for high-quality Keck 
spectra, the first is probably negli- 
gible, while the latter two are com- 
parable. 

How do the thermal broaden- 
ing scale and the baryon-smoothing- 
scale change with reionization his- 
tory? For a fixed redshift of obser- 
vation (say 2 = 3), as one raises the 
redshift of reionization, the temper- 
ature at z = 3 becomes lower (Hui 
k, Gnedin 1997)| and the thermal 
broadening scale becomes smaller. 
The effect on the baryon-smoothing- 
scale is more subtle. It is true that 
the Jeans scale, like the thermal 
broadening scale, is proportional to 
VT, and so lowering the tempera- 
ture lowers both. However, as shown 
by Gnedin & Hui (1998), the true 
baryon-smoothing-scale is in fact not 
given by the Jeans scale, but gener- 
ally given by something smaller. In 
fact, right before reionization occurs, 
the baryon-smoothing-scale is very 
small, and afterwards, the baryon- 
smoothing-scale does not suddenly 
jump up to the conventional Jeans 

2 T is quite independent of Jjjj as 
long as H is sufficiently ionized, but 
does depend somewhat on JhcII- 



scale, but only slowly catches up with 
it. This means that raising the red- 
shift of reionization could have the 
opposite effect of allowing a larger 
baryon-smoothing-scale by allowing 
more time for it to catch up with the 
Jeans scale. Therefore, exactly how 
raising the redshift of reionization af- 
fects the b-distribution requires a de- 
tailed calculation. 

This is a subject of much current 
interest, particularly because of the 
work of Bryan et al. (1998) (see also 
Haehnelt & Steinmetz 1998, Theuns 
et al. 1998) who, after carefully in- 
vestigating the effect of numerical 
resolution and box-size on the b- 
distribution, found that the canon- 
ical as = 0.7 SCDM model predicts 
a b-distribution that has too many 
narrow lines than is observed, if the 
universe reionizes by z ~ 6. The in- 
teresting questions are: I. changing 
the fluctuation amplitude will likely 
shift the b-distribution, but which 
way will it go? - the arguments lead- 
ing to eq. (^) indicate that lower- 
ing the amplitude would move the 
b's up (because lowering the am- 
plitude means more suppression of 
sharp fluctuations or narrow lines), 
but as pointed out by Hui & Rut- 
ledge (1998), nonlinear corrections 
could reverse this trend; II. reion- 
ization history will no doubt affect 
the b-distribution, but as pointed out 
above, which direction it will go is 
not obvious and requires a detailed 
calculation; III. as noted by Hui & 
Gnedin (1997), the mean tempera- 
ture of the intergalactic medium in- 
creases with Cj, (see also Bryan et 
al. 1998) and decreases with f2 m ; it 
would be interesting to see whether 
changing these would fix the discrep- 
ancy of Bryan et al.. 

It is clear from the above discus- 
sion that the inference on reioniza- 
tion history from the b-distribution 
will depend on assumptions made 
about the cosmological density pa- 
rameters and the power spectrum 
normalization. To isolate the effect 
of reionization history from the ef- 
fect of power spectrum normaliza- 
tion, one possibility is to consider 
the lowest b's, which according to 
our picture, should correspond to 
the minimum size imposed by ei- 
ther the thermal-broadening scale, 
or the baryon-smoothing-scale. How- 
ever, one should keep in mind the 



4 



Lam Hui 



possible complication that not all fit- 
ted absorption-lines arise from peaks 
in r, e.g. some narrow lines might be 
introduced to "fill in" the wings of 
absorption peaks which do not nec- 
essarily have Voigt-profile shapes, in 
which case the physical meaning of 
the widths of these lines is unclear. 
One should check for this possibil- 
ity using simulations, or even develop 
other characterizations of the line- 
width which are less prone to system- 
atics of this sort. 

3.2 The Primordial Mass Power 
Spectrum 

This is an area pioneered by Croft 
et al. (1998), who showed that the 
linear mass power spectrum can be 
reliably recovered from the power 
spectrum of the transmission. The 
idea works as follows. The trans- 
mission power spectrum Pf along 
the line of sight and the three- 
dimensional mass power spectrum 
P p are related to each other on large 
scales by (Hui 1998): 




BW{k\\,k)P p {k)kdk(A) 



where P^ is the fourier trans- 
form of the two-point correlation 
(5f(0)S f (u)) (S f = (/-/)// where 
/ is the mean transmission), fcy is 
the velocity-wave-vector along the 
line of sight and k is the mag- 
nitude of the corresponding three- 
dimensional wave-vector. The ker- 
nel W(kpk) describes the redshift- 
distortion of the power spectrum. B 
is a constant which is determined 
by the nonlinear transformation from 
density to the transmission exp(— r) 
(eq. pj). Fortunately, the only free- 
parameter that enters into the de- 
termination of the constant B is the 
constant A in eq. (jj|) which can 
be fixed by matching e.g. the ob- 
served mean transmission (Croft et 
al. 1998). It is important to real- 
ize that we have made use of the 
facts that a has a small range (eq. 
|lj; jj^), that the uncertain thermal- 
broadening and baryon-smoothing- 
scales only affect the transmission 
correlation on small scales, and that 
the dark matter and baryon distribu- 
tions trace each other on large scales. 



It is expected eq. (^) holds on large 
scales with P p being the linear power 
spectrum. Hence, once the distortion 
kernel is fixed (which is predicted 
by gravitational instability, with de- 
pendence on the cosmological den- 
sity parameters H's), the primordial 
mass power spectrum can be recov- 
ered from pi (k) by inverting an es- 
sentially triangular matrix propor- 
tional to the distortion kernel (Hui 
1998). 

There are two main advantages of 
this way of measuring the linear mass 
power spectrum. First, unlike in the 
case of galaxies, the "biasing" rela- 
tion between density and the observ- 
able (the transmission) is known ex- 
actly here (aside from the parame- 
ter A which can be fixed using inde- 
pendent observations). Second, this 
provides a direct probe of the linear 
mass power spectrum on small scales 
(fc<0.02s/km ~ 4/i v / n^7Mpc- 1 , for 
z ~ 3; see Croft et al. 1998), which 
are out of the reach of galaxy surveys 
at low redshifts (unless one corrects 
for the nonlinear evolution). This is 
simply because the nonlinear scale 
becomes smaller at higher redshifts. 

Exciting (and beautiful!) first re- 
sults of the application of the above 
techniques to the observed forest 
were reported recently by Croft et al. 
(1998b), Weinberg et al. (1998) (see 
also contribution to this volume by 
Weinberg). Here, let us list a few is- 
sues that deserve some thought and 
perhaps further investigation. I. The 
distortion kernel W above can be 
predicted using linear theory, but it 
is quite possible that the highly non- 
linear transformation from the den- 
sity to the transmission will alter its 
behavior, even at very large scales 
(McDonald & Miralda-Escude 1998, 
Hui 1998). This should be carefully 
checked using simulations. II. There 
is an upper limit to the scale above 
which we cannot reliably recover the 
mass power spectrum. It is set by the 
continuum, which is known to have 
long range fluctuations, and which 
can only be estimated up to a limited 
accuracy. This limit is around k ~ 
0.002s/km, but should be checked 
carefully using high resolution ob- 
servations. III. The whole inversion 
procedure discussed above relies on 
the fact that eq. (|lj) holds for most 
of any quasar spectrum, aside from 
regions of strong absorption such as 
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those caused by collapsed halos. It 
would be useful to have an estimate 
of how such regions affect the recov- 
ery procedure. IV. Numerical sim- 
ulations are typically used to effec- 
tively fix the constant B in eq. (^). 
It will be useful to have a study 
of the minimal box-size and reso- 
lution required for the problem at 
hand. For instance, while one expects 
that resolution can be sacrificed as 
long as one is interested in large 
scale fluctuations, resolution does af- 
fect the fixing of the amplitude A 
using the observed mean transmis- 
sion (Croft, priv. comm.). Analytical 
methods to fix B would be very use- 
ful. V. The temperature-density re- 
lation mentioned in is expected to 
have a scatter. It would be good to 
have an idea by how much it could 
modify the transmission power spec- 
trum, or in other words, the "bi- 
asing" relation between density and 
transmission. 

3.3 The Cosmological Energy 
Contents 

Here, we discuss a version of a 
test proposed by Alcock & Paczyhski 
(1979; AP hereafter), which is par- 
ticularly sensitive to the presence 
of the cosmological constant A, or 
more generally, a component of the 
cosmological energy contents which 
has negative pressure, let us call it 
Q, with an equation of state p = 
wp (w < 0). The ideas presented 
here have been considered by several 
groups recently (Croft 1998, Seljak 
1998, McDonald & Miralda-Escude 
1998, Hui, Stebbins & Buries 1998). 
AP observed that an object placed at 
a cosmological distance would have 
a definite relationship between its 
angular and redshift extents, which 
is cosmology-dependent. Consider an 
object with mean redshift z, and an- 
gular size 9. Its transverse extent in 
velocity units is 

u ± (6) = -Z-D A {z)0 . (5) 

1 + 2 

Here H is the Hubble parameter at 
redshift z, and D A (z) is the an- 
gular diameter distance (Weinberg 
1972). For spherical objects the ra- 
dial and transverse extents are equal, 
but more generally if the object is 
squashed radially by a factor a s , 



the radial extent is Wii = fn - ^ = 

II 1 + z 

a s u±. Here c is the speed of light 
and uj_, U|| < c is assumed. This 
parameter Uj_/6 depends on all the 
different cosmological density pa- 
rameters f2 m (matter), (curva- 
ture) and Qq (Q), but is particu- 
larly sensitive to Hq. It is signifi- 
cantly lower for a Q-dominated uni- 
verse (for w$ — 1/3) than for a 
no-Q-universe. (see Fig. 1 of Hui et 
al. 1998) Qualitatively, one can un- 
derstand this as follows. In a Q- 
dominated universe, it is well known 
that a significantly larger radial- 
comoving-volumc is associated with 
a given redshift range Az, compared 
to a no-Q-universe i.e. cAz/H(z) is 
larger for a Q-dominated universe. 
Hence, the angular-extent of a spher- 
ical object of a given redshift-extent 
Az, would appear to be larger for 
a Q-universe: 8 = cAz/H(z)/D A (z) 
would be larger i.e. a smaller uj_(8). 
It turns out for an open universe with 
no Q, H(z) and Da{z) roughly bal- 
ance each other to make uj_ (8) simi- 
lar to that of a fiat-no-Q-universe. 

In the case of the Lya forest, 
the "object" to use is the two-point 
correlation function, whose "shape" 
(ct B ) is not spherical because of 
redshift-anisotropy induced by pe- 
culiar motion. We cannot observe 
the full three-dimensional correla- 
tion directly; instead we can mea- 
sure the one-dimensional correlation 
along a line of sight, and the cross- 
correlation between two close-by 
lines of sight, or their Fourier coun- 
terparts: the auto- and the cross- 
spectra. A comparison of the two 
provides a new version of the AP 
test. More precisely, given the ob- 
served transmission Dower spectrum, 
one can invert eq. <M) to obtain the 
mass power spectrum, and then use 
the following equation to predict the 
cross-spectrum P^, between two lines 

of sight with separation 8 (P^ is 
the fourier transform of the cross- 
correlation between the two lines of 
sight a and b: {5 a f (0)S b . («))): 

P OO 

P f x (k h 8)= I BW(k }] /k,k) (6) 

P p (k)J [k ± u ± (8)]kdk 

where Jo is the spherical Bessel func- 
tion. This prediction is cosmology 
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dependent mainly through the pa- 
rameter and also through the 
distortion kernel W(ku/k,k). The 
latter generally depends on an ef- 
fective bias parameter (Hui 1998), 
aside from the cosmological density 
parameters of interests, which can 
fortunately be determined from sim- 
ulations because the exact "biasing 
relation between density and the ob- 
servable (the transmission) is known 
(eq. JlJ), unlike in the case of galax- 
ies (where the AP test has been con- 
templated by other authors). A com- 
parison of the predicted with the 
observed cross-spectra, for a given 
observed auto-spectrum, provides a 
version of the AP test: assuming 
the wrong cosmology will result in 
a wrong prediction for the cross- 
spectrum. Note that in this test, 
there is no need to fix the constant B 
because it enters into the auto- and 
cross-spectra in the same way. 

According to Hui et al. (1998), 
to reach a 4 — a level discrimina- 
tion between e.g. the Q m = 0.3 - 
Hj. = 0.7 universe and the Q m = 0.3 
- f^A = 0.7 universe, only 25 pairs 
of quasar spectra, at angular sepa- 
rations 0.5' — 2', would be required. 
There are roughly 10 such pairs of 
quasars with existing spectra at the 
above angular separations, or slightly 
larger, and at zijl (see e.g. Crotts & 
Fang 1998 & ref. therein). Upcoming 
surveys such as the AAT 2dF and 
SDSS are expected to increase this 
number by at least an order of mag- 
nitude, and to higher redshifts. 



4 Tests and More 

The cosmological utility of the 
(low column density) forest relies 
very much on the simple relation be- 
tween the optical depth or transmis- 
sion and the density+velocity fields, 
which arises from the fact that I. 
the distribution of the latter is de- 
termined by gravitational instability 
alone on large scales, and by baryon- 
smoothing on small scales (e.g. ex- 
plosions do not play an important 
dynamical role); II. the ionizing 
background J does not have signif- 
icant spatial fl uctu ations on scales of 
interest (see §3.2). Here, we loosely 



rent simulation-inspired work on the 
forest, which makes use of the above 
assumptions, to explain the observed 
b and column-density distributions 
should be counted as a vindication 
of the paradigm. After all, this body 
of work is based on structure for- 
mation models constructed to match 
observations other than those of the 
forest. Moreover, there are theoreti- 
cal reasons to believe that J fluctu- 
ations should be small on the scales 
we are interested in (e.g. Croft et al. 
1998b). Nonetheless, it is important 
that we find alternative ways to test 
our assumptions observationally, as 
we continue to look for new applica- 
tions of the forest in cosmology. 

Here, we discuss three possible 
tests that involve only the forest, and 
one cross-test with another high red- 
shift observation. 

I. Use double (or multiple) lines 
of sight. There are in fact two dif- 
ferent tests here. First, for quasar 
pairs that are very close together e.g. 
lensed pairs, one can check whether 
the cross-correlation between the for- 
est in the two lines of sight is as 
strong as one would expect based 
on the smooth-fluctuation-paradigm. 
More specifically, for lensed quasar 
pairs which are typically of the or- 
der of arc-second separation, the dis- 
tance between the two lines of sight 
is sub-kpc, which is much smaller 
than the baryon-smoothing-scale ex- 
pected for a gas of T ~ 10 K. 
This means one expects 100% cor- 
relation between these two lines of 
sight through the forest. Exactly 
such a behavior has been observed 
by Rauch (1997). This shows there 
is no fluctuation in the ionizing back- 
ground, nor is there complicated mo- 
tion in the forest due to explosions, 
at least on these small scales. For 
quasar-spectra at larger separations, 
on the other hand, assuming a given 
set of f2 's. o ne can turn the proce- 
dure in 1 3.3 around, and test for the 



refer to this whole set of assump- 
tions as the smooth-fluctuation- 
paradigm. The ability of most cur- 



redshift-anisotropy due to peculiar 
motion. This is a robust prediction of 
gravitational instability. Any other 
source of fluctuations, such as that 
due to a fluctuating background, will 
act as additional sources of "bias- 
ing" , and change the prediction for 
the redshift-anisotropy. 

II. Use higher order statistics. It 
is well known that gravitational in- 
stability gives robust predictions for 
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higher order statistics such as the 
skewness (5 3 }/(5 2 } (Peebles 1980), 
where 8 is the overdensity in mass. 
Just as in the case of galaxies, 
where measuring the skewness for 
the galaxy distribution ((<5^) / {<5g} 
where 8 g is the galaxy overdensity) 
provides information on the biasing 
relation between mass and galaxy 
distributions, measuring the skew- 
ness of the transmission should pro- 
vide a test of the "biasing" relation 
between transmission and mass den- 
sity (eq. jj|). Any additional compli- 
cation due to say a fluctuating ion- 
izing background is going to change 
that relation and hence the predic- 
tion for the higher moments. Ef- 
forts are under way to calculate these 
quantities. As far as measurements 
are concerned, it is important to 
guard against possible estimation- 
biases (Hui & Gaztanaga 1998). It 
might also well be the case that 
other higher order measures are bet- 
ter suited for the purpose at hand, 
because of the particular highly non- 
linear transformation from density to 
transmission we have here (especially 
the exponentiation: exp(— r)). 

III. If supernova explosions play 
an important role in the distribution 
of the gas that makes up the for- 
est, we should be able to see their 
remnants: metals. Clever efforts to 
detect metals in Afjji < 10 14 cm — 2 
forest have so far only yielded up- 
per limits (see Lu et al. 1998). How- 
ever, a carbon-to-hydrogen ratio of 
the order of 10 — 2 5 has been detected 
in systems with Nm ~ 10 14 ' 5 cm — 2 . 
It is unclear, though, what the cor- 
responding volume-filling factor of 
these slightly-metal-enriched regions 
is. This aspect of the forest certainly 
deserves closer scrutiny. ^| 

Finally, it is appropriate at this 
conference to make some connec- 
tion with another probe of the 
high-redshift universe: the Lyman- 
break-galaxy survey (Steidel et al. 

3 Note added after conference: 
Cowie & Songaila (1998) recently re- 
ported results that are in apparent 
contradiction with Lu et al. (1998). If 
the results hold up, this implies some 
amount of star formation occured at 
early times, but has stopped by z ~ 3 
in order not to violate point I above. 
Its impact on the forest remains to 
be investigated. 



1998). As noted by several authors 
(Nusser & Haehnelt 1998, McDon- 
ald & Miralda-Escude 1998, Croft et 
al. 1998), information obtained from 
the Lyman-alpha forest can be fruit- 
fully combined with observations of 
the Lyman-break objects, which are 
observed at similar redshifts (z ~ 3). 
For instance, once the mass power 
spectrum is inferred from the forest, 
a comparison with the clustering of 
the Lyman-break objects would im- 
mediately yield a measure of the bias 
of the Lyman-break galaxies. On the 
other hand, the bias of the Lyman- 
break galaxies can be obtained from 
observations of the Lyman-break ob- 
jects themselves: either by a direct 
measurement of their masses (us- 
ing high resolution infrared observa- 
tions, which should be feasible in the 
near future), or by measuring the 
skewness of the galaxy-distribution. 
Note that to deduce from either 
quantity the bias (as defined by the 
square root of the ratio of the galaxy 
power spectrum to mass power spec- 
trum), one needs a model for how 
these galaxies form i.e. how they are 
formed in relation to dark matter ha- 
los (Mo et al. 1998). 

With the above information in 
hand, we have two independent mea- 
sures of the mass power spectrum at 
z ~ 3, one from the Lyct forest, the 
other from the Lyman-break-objects 
making use of the deduced bias pa- 
rameter. They probe the clustering 
at different scales, and so one can use 
them to constrain the shape of the 
power spectrum; or, if one assumes 
a particular cosmological model with 
a definite power spectrum shape, 
one can test for consistency. More- 
over, the latter approach also al- 
lows two independent measurements 
of f! m (with minor dependence on 
the other parameters), by compari- 
son with the cluster-normalized mass 
power spectrum today (see e.g. Gi- 
avalisco et al. 1998, Adelberger et al. 
1998, Weinberg 1998). Investigations 
along these lines are being pursued. 
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